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We propose a method to evaluate the entropy density and entropy flux in a vacuum gap between 
two half-spaces that takes into account influence of near-field effects, i.e., interference, diffraction, 
and tunneling of waves. The method developed is used to determine the maximum work that can 
be extracted through near- field radiative transfer between two half-spaces at different temperatures. 

PACS numbers: 44.40.+a,42.25.-p,05.70.Ln,05.07.-a 



The inception of modern quantum physics can be 
traced back to Planck's pioneering work on blackbody 
radiation [1 . Central to Planck's work is the thermo- 
dynamic analysis of thermal radiation in a cavity, which 
requires knowledge of energy, momentum, and entropy of 
photons. Planck's analysis is restricted to the case when 
near-field effects (NFE), i.e., the collective influence of 
diffraction, interference, and tunneling of waves, are ab- 
sent m . Though we know much about energy (leading to 
radiative heat transfer) and momentum (leading to van 
der Waals and Casimir forces) associated with near field 
radiative transfer (NFRT), little, if any, is known about 
entropy [2 . The purpose of this letter is to determine en- 
tropy density and entropy flux in a planar vacuum cavity 
between two half-spaces at different temperatures when 
the separation between them is small enough that NFE 
are important. Using the entropy flux, the maximum 
work that can be extracted through NFRT between two 
half-spaces at different temperatures is determined. 

Planck showed that the far-field radiation entropy in- 
tensity is given by [1 : 

(^) = [(1 + M^) ln(l + M^) - In M^] (1) 

where /c^ is the Boltzmann constant, uj is the angular 
frequency, c is the speed of light in free space, = 

M(cj, Q) = ^^^^] ^\ o , 27rfi is the Planck constant, and 

I(uj^ (l) is the spectral radiation intensity of a ray in the 
direction defined by unit vector Vt. /3 = 1 for polarized 
radiation and /3 = 2 for unpolarized radiation. The spec- 
tral entropy density is given by = (1/c) J l^{Cl)dQ. 
von Laue [3j and other researchers [4 investigated the 
effect of partial coherence on the entropy of radiation. 
Petela and others O |6] used the concept of availability or 
exergy to determine the maximum work that can be ex- 
tracted from solar radiation. All of the above mentioned 
works [3]-[6] rely on Eq. [l]for entropy intensity. Under- 
lying the derivation of Eq. [l] is the assumption that the 
electromagnetic local density of states (LDOS) and veloc- 
ity of radiation are independent of position and given by 
their values in free space. The implications of the break- 
down of the above-mentioned assumptions when NFE are 



present, to the best of our knowledge, have not been in- 
vestigated. 

Consider a planar vacuum cavity of thickness Iq = 
\zr — zl\ between two parallel half-spaces L and as 
shown in Fig. [l] at temperatures Tl and Tr respec- 
tively. Each half-space is assumed to be isothermal. Half- 
spaces with planar thin films do not introduce any con- 
ceptual difficulties and are taken into account via gen- 
eralized Fresnel reflection coefficients Rql and Ron- Be- 
cause of translational symmetry in the xy plane, all quan- 
tities pertinent to this letter only depend on z. The 
influence of NFE is captured by modeling thermal ra- 
diation using Rytov's theory of fluctuational electrody- 
namics (FE) [7 . In FE, the cross power spectral densi- 
ties of the fluctuating current sources are related to the 
electromagnetic properties and temperature through the 
fluctuation-dissipation theorem [8 . The ensemble aver- 
aged Poynting vector and energy density can be writ- 
ten in terms of the dyadic Green's functions (DGFs) of 
the vector Helmholtz equation, which are well known for 
layered media [9 . Though fluctuational electrodynam- 
ics provides us with a framework for determining energy 
density and Poynting vector, it does not tell us how to 
compute entropy. 

Entropy can be calculated by determining the number 
of microscopic states into which a given number of pho- 
tons can be distributed. Consider a thin film of area A, 
which can be made as large as necessary, lying between 
coordinates z and z-\-dz in the vacuum gap. We first con- 
sider the case when Tl 7^ Tr. The entropy of electromag- 
netic waves emitted from half-space L with microscopic 
states in the range /i to ji ^ dji within this thin film is 




FIG. 1. (Color online) Two half-spaces L and R separated by 
a vacuum gap of thickness /o- 
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given by sl{ii^ z)diiAdz and obtained from the relation: 
sl{ii^ z)diiAdz = /cb In SWl- Here fi is the set of variables 
defining the space containing the microscopic states, 
which is unidentified as yet, and 5Wl is the number of 
ways in which riLif^^ z)djiiAdz photons can be distributed 
across the number of accessible microscopic states, rep- 
resented by pa{lJi^z)diiAdz. nL{ii^z) and pa{/J^,z) are the 
number density of photons and accessible microscopic 
states per unit volume in the /i space. The total num- 
ber of available microscopic states is related to the LDOS 
p{uj, z) and is given by z)diiAdz^ of which only a frac- 
tion pa{p^ z)/ p{p^ z) is accessible to the photons from L. 

OWl is given Dy [iU ^nL{^Ji,z)d^lAdz)\{pa[^Ji,z)d^JLAdz-l)\' ^ 

similar expression exists for SWr. Taking the two inde- 
pendent polarizations into account, we can write (assum- 
ing that pa{p^ z)dpAdz 1): 



s\^\p,z) =kBp'^Kp.z) 



\p,z) 



In 



{p,z) 



Pa\p,z) 



1 



(2) 



where h = L^R refers to contributions from half-spaces L 
and R respectively, and j = s^p refers to 5 or p polariza- 
tion. Since the electromagnetic waves emitted from the 
two half-spaces are incoherent, the total entropy density 
at any location is given by s{z) = J dp ^ {s^^\p^z) + 



{p,z)) [ir 



Aj) 



can be determined from its relation to 
(z), the energy density in the vacuum gap of s or 
p polarized waves from half-space h. The method used 
by Antezza et. al. can be used to determine u^^\z) 
[13]. n^^\p^z) and u^^\z) are related as follows: 



(^) 



dphujn^l\p^ z)-\- / dphjbJirSi\p^ z) (3) 



IPW JEW 

where PW{EW) refers to propagating (evanescent) 
waves. The integrals Jp^ and will be defined 



shortly. n^i\p^z) is given by (Eq. 
4b] for EW): 



4a 



for PW and Eq. 
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where kp is the magnitude of the in-plane wave vector 

= kpl<.p = kxX + kyy^ ko = cj/c, /c^ + k'^Q 
{kp < ko), and kp 



kl for PW 



Pio = k'o ioiEW {kp>ko). 
Constant ko surfaces, which are spheres for PW {k^ + 
kl = k^) and hyperboloids for EW {k^ - fS^ = k^), 
are shown in Fig. [2|l and Fig. respectively (only 
kz^ Pz > ^ portion is shown). Since the expression for en- 
ergy density (or stress tensor) is an integral over the in- 
plane wave vectors and frequency, we can identify p as the 
space spanned by {kp G (0, oo),0 G (0,27r),/co G (0,oo)} 
{ko is used instead of cj). Consider a differential patch 
on a constant energy surface (blue patch in Fig. [2|i and 
Fig. ^p) whose projected area in the kx — ky plane is 
dkxdky or kpdkpdcj). The area of this blue patch is given 

by dS = - — ^ where Oz{p) is the angle between the 

I cos Oz{p)\ 

surface normal to the constant energy surface and the kz 
or axis, cos 6 z{p) = kzo/ko for PW and cos 6 z{p) = 



Pzo/ y PzO ~^ kp EW. dp is an infinitesimal volume 
element of area dS on a constant ko surface and thick- 
ness dko perpendicular to it. The magnitude of dp is 
dkodS. Jpy^ dp and J^^^ dp represent triple integrals 
over the domains {(j) G (0,27r),A:p G (0, /co),/co G (0,oo)} 
and {(j) G (0, 27r), kp G (/cq, oo), ko G (0, oo)} respectively. 

Since pa\p^z) is related to p^^\p,z), the integral 
expression for p{uo, z) will be used to first determine 

p^^\p,z). In vacuum p(a;,r) = :^^Tr(Ge(r, r; cj) + 
Gm(r, r;*^)) [14 , where Tr(G) is the trace of G, and 
Ge(r, r; cj), Gm(r, r; cj) are the electric and magnetic 
DGFs. p^^\uj,z)duj [15^ is given by: 



p^^\uj, z)d(jj 



dpp^^\p,z) 



PW 



dpp^^\p,z) (5) 



EW 



where j'p^ dp dko Jpy^ dS and dp dko Jew 
are double integrals over the domains {(j) G (0,27r),/Cp G 
(0, ko)} and {(j) G (0 , 27r), kp G (/cq, oq)| respectively, and 



p(^)(/i,z) is (Eq. l6aJfor PW and Eq. l6b]for EW): 
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From the dispersion relation, we see that kzo = 
± y^/^o — and /^^o = ^yj^p — k^. Only the posi- 
tive sign for kzo and /^^o is valid in the expressions for 
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n^i\li^z) (Eq. In contrast, both signs of the square 
roots are vahd in Eq. |6|for p^^\iJi^ z). Using the property 
that RoL{-kzo) = ^ol(^2o) and RoR{-kzo) = RQ^ik^o) 
[9l see Eq. 2.7.3], it can be seen that the same value 
of p^^\ii^ z)dii is obtained irrespective of the choice of 
sign for the square root, i.e., /c^,/^^ > and A:^,/^^ < 
portions of the ji space, corresponding to states accessi- 
ble to waves originating from L and R respectively, con- 
tribute equally to p^^\p^ z)dp. Hence, for thermal non- 
equilibrium condition, pa\p^z) = p^^\p^ z)/2. The case 
of thermal equilibrium is discussed later. 

The polarization dependent energy flux in the z direc- 
tion due to sources in L, E^^ (it is independent of z), 
can be written as 



d/j, cos 0^ (m)©(w, Tl)K^^^ iij) 



(7) 



where 6(w,Ti,) = hw / {exp{fkj / kBTi) — 1), and 
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Biehs et. al. [16 and Ben-Abdallah and Joulain [17] have 
interpreted {^ii^ / c)K^^\p) as a generalized transmissiv- 
ity between the two half-spaces L and R. Here, we in- 
terpret {uj'^/c^)Q{uj,Tl)K^^\p) as a generalized spectral 
radiation intensity^ valid for PW as well as EW . The 
concept of radiation intensity (/) in classical theory of 
thermal radiation is associated with the power contained 
in a cone of dVt. For PW, dS/k^ can be associated with 
the solid angle in the direction of propagation because 
the constant /cq surface happens to be spherical. The 
solid angle interpretation of intensity is invalid for EW. 
K^^\p) on a constant k^ surface for PW and EW (the 
entire surface is not shown) are shown in Fig. |2]3 and Fig. 
[2]i respectively (simulation performed for s polarization, 
/o = 5 /im, 27r//co = 10 /im, = sr = 2.2 -|- O.Olz). In 
each of the two figures, the length of the arrows are pro- 
portional to the magnitude of K^^\p). The arrows are 
perpendicular to the constant k^ surface. As /cp ^ oo, 
the magnitude of the arrows decreases exponentially, in 
accordance with Eq. IS] 
Xi) 



M can be also expressed as 



E 



dphwn'l^ (/i, z)v[^l (/i, z) 



(9) 



where v[^l{p^z) is the z component of the polarization 
dependent local velocity of energy transmission associ- 
ated with photons from L. Since K^i^n) is the same for 
both L and R half-spaces, and cos6z{p) is opposite in 
sign, E^^ is opposite in direction to E^^ . At equilibrium, 
they cancel each other to yield zero net radiative trans- 
fer between the two half-spaces. Hence, at equilibrium 




FIG. 2. (Color online) Constant frequency surfaces for (a) 
PW, and (b) EW. The surfaces are symmetric about the kx — 
ky plane. In (c) and (d), the arrows are proportional to and 
point in the direction of the generalized radiation intensity 
for (c) PW, and (d) EW. 



be shown to be: 
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fc2 2U(R: 



by comparing Eq. |4] and Eq. [Sj We have confirmed nu- 
merically that, as one would expect, v[^l{p, z) < c for all 
p. The entropy fiux, Sl{z), associated with this energy 
fiux is given by 



Sl{z) 



dp 



So) 



(11) 



Eji{z) and Sr{z) can be obtained by the same procedure 
if Tl 7^ Tr. Unlike energy fiux within the cavity which is 
independent of z, though any two waves with different p 
are incoherent, interference effects due to multiple refiec- 
tions of the same wave lead to a z dependence of entropy 
fiux [H. 

The maximum work that can be extracted from NFRT 
between the two half-spaces (rejecting radiation at sink 
temperature Tr) in Fig. [l] is given by Wmax = {El - 
TrSl) — {Er—TrSr) |6|. In this expression, both and 
Sr are evaluated at z = zr. The maximum efficiency of 
work extraction is then given by rjmax — Wmax /El- Park 
et. al. [18 analyzed the performance of a thermopho- 
tovoltaic (TPV) converter using W emitter at 2000 K 
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a generalized spectral radiation intensity which can be 
used even when NFE are present. Though the concepts 
developed here are used to analyze energy conversion us- 
ing NFRT, it can also be used to better understand the 
thermodynamics of surface wave-based laser cooling [2T] 
and thermal non-equilibrium Casimir interactions [22] . 

This work is funded partially by National Science 
Foundation Grant CBET-0853723. 



FIG. 3. (Color online) Energy and entropy flux (left y axis, 
hollow markers) and conversion efficiency (right y axis, solid 
markers) as a function of vacuum gap for W-Ino.i8Gao.82Sb, 
compared to results in Ref. 18 . 

and Ino.i8Gao.82Sb photovoltaic cell at 300 K. Using op- 
tical data for W J9] and Ino.i8Gao.82Sb [20 , we compute 
El-, Sl, Er^ and Sr, of which El and TrSl are shown in 
Fig. |3l when L is W at 2000 K and R is Ino.i8Gao.82Sb 
at 300 K. Both energy flux and entropy flux increase 
with decreasing gap, characteristic of tunneling due to 
evanescent waves. We also plot rjmax as a function of 
gap and compare with the efficiency predicted by Park 
et. al. [18]. As observed by Landsberg and Tonge [6], 
though the high values of r]max are usually unattainable, 
the utility of r]max is to impose an upper limit on efficien- 
cies of all models for energy conversion, including TPV 
conversion, involving the same materials, configurations, 
and temperatures. 

The case of Tl = Tr = T has to be treated differ- 
ently since emission from both half-spaces are at the same 
temperature. When Tl = Tr^ electromagnetic waves in 
the vacuum cavity cannot be distinguished as originat- 
ing from L OT R because their temperatures are equal. 
Hence, in thermal equilibrium, it is n^i\/J^^ z) -\-n^^\/j.^ z) 
rather than n^i\/j.^z) or n^^\fi^z) that determines the 
entropy. At the same time, pa\l~i^z) = p^^\fi^z) since 
kz^Pz ^ and kz^Pz ^ portions of the dispersion 

surface are accessible now. Since ^~^'^..w — ^ = 

0(0;, T) at thermal equilibrium, the contribution to 
entropy density is given by p*^-^^ (/i, z)7(a;, T), where 

li'-^T) = fcB[ e.p(tg)-i - - exp(-^)]] [T. 

Of this, a fraction Q^^T)p^'j)'{il z) froi^i half-space L 
and the remainder from R. The equilibrium entropy 
flux in the cavity due to thermal sources within L is 

given by S^j^ = J dfi ^^l\^^ ^)^i^L{^^ ^)q0X) 

{S^ is independent of z). Since n^^\/jj^ z)v^J^{fi^ z) = 
—n^i\/i^ z)v^^l{li^ z) at thermal equilibrium, S^j^ and S^^ 
cancel each other and the net entropy flux is zero. 

In summary, we elucidate a method to determine en- 
tropy density and entropy flux between two half-spaces 
when near field effects are important. We also identified 
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